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Abstract. This paper studies the Levy model of the optimal multiple-stopping problem arising in the 
context of the timing option to withdraw from a project in stages. The profits are driven by a general 
exponential spectrally negative Levy process. This allows the model to incorporate sudden declines of the 
project values, generalizing greatly the classical geometric Brownian motion model. We solve the one- 
stage case as well as the extension to the multiple-stage case. The optimal stopping times are of threshold- 
type and the value function admits an expression in terms of the scale function. A series of numerical 
experiments are conducted to verify the optimality and to evaluate the efficiency of the algorithm. 
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Consider a firm facing a decision of when to abandon or contract a project so as to maximize the total expected 
future cash flows. This problem is often referred to as the abandonment option or the contraction option. A typical 
formulation reduces to a standard optimal stopping problem, where the uncertainty of the future cash flow is driven 
by a stochastic process and the objective is to find a stopping time that maximizes the total expected cash flows 
realized until then. A more realistic extension is its multiple-stage version where the firm can withdraw from a 
project in stages. 

In a standard formulation, given a discount rate r > and X t = x + (p — \o~ 2 )t + o~Wt for a standard Brownian 
motion W, p € K and a > , one wants to obtain a stopping time r of X that maximizes the expectation 



The profit collected continuously is modeled as the geometric Brownian motion e ' less the constant operating 
expense 5 > 0. The value corresponds to the lump-sum benefits attained (or the costs incurred) at the time 

of abandonment. Here a technical assumption r > p is commonly imposed so that the expectation is finite and 
the problem is non-trivial. The problem is rather simple mathematically; it reduces to the well-known perpetual 
American option (or the McKean optimal stopping problem). An explicit solution can be attained even when X is 
generalized to a Levy process (see, e.g., Mordecki [26]). 



1. Introduction 



(1.1) 
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In this paper, we generalize the classical model by extending from Brownian motion to a general Levy process 
with negative jumps (spectrally negative Levy process), and consider the optimal stopping problem of the form: 



(1.2) supE 



[ T e- rt f(X t )dt + e- rT g(X T )l {T<oo} 
Jo 



We obtain the optimal stopping time as well as the value function for the case / is increasing and g admits the 
form g(x) = K — Yl%=i cie a%x for some strictly positive constants c and a. We also show the optimality among 
all stopping times of threshold type (see (2.4) below) when g is relaxed to be a general decreasing and concave 
function. 

We further extend to the multiple-stage case where one wants to obtain a set of stopping times {r( m ); 1 < m < 
M} such that = < rW < • • • < r( M ) a.s. and achieve 



M 



,. T ( m ) 

(1.3) sup V E 

r (l)<...< T (M) ^ 



/ e- rt F m (X t )dt + e-^ (m) ffm (X r(m) )l {r(m)<oo} 

J T (m-l) x s 



when g m and f m := F m — F m+ i (with Fm+i = 0), for each 1 < m < M, satisfy the same assumptions as in the 
one-stage case. The multiple-stopping problem arises frequently in real options (see e.g. [14]) and is well-studied 
particularly for the case X is driven by Brownian motion. In mathematical finance, similar problems are dealt in 
the valuation of swing options [10, 11]. 

Although the use of the geometric Brownian motion is fairly common in real options, empirical evidence sug- 
gests that the real world is not Gaussian, but with significant skewness and kurtosis (see, e.g., [8, 13, 30]). Dixit 
and Pindyck [14] considered the case with jumps of a fixed size with Poisson arrivals. Boyarchenko and Leven- 
dorskii [9] considered the EPV approach for a general Levy process satisfying the (ACP)-condition (with a focus 
on exponential-type jumps for illustration); they solved a related multiple-stage problem with g being constant. The 
Levy model is in general less intractable than the continuous diffusion counterpart, especially when the lump-sum 
reward function g is not a constant. When jumps are involved, the process can potentially jump over a threshold 
level, requiring one to compute the overshoot/undershoot distributions that depend significantly on the form of the 
Levy measure. Technical details are further required when it has jumps of infinite activity or infinite variation. 

In this paper, we take advantage of the recent advances in the theory of the spectrally negative Levy process 
(see, e.g., [7, 21]). In particular, we use the results by Egami and Yamazaki [17], where we obtained and showed 
the equivalence of the continuous/smooth fit condition and the first-order condition in a general optimal stopping 
problem. Unlike the two-sided jump case, the identification of the candidate optimal stopping time can be con- 
ducted efficiently without intricate computation. The resulting value function can be written in terms of the scale 
function, which also can be computed efficiently by using, e.g., [18, 29]. The extension to the multiple-stage can 
be carried out without losing generality. The resulting optimal stopping times are of threshold type with possibly 
simultaneous jumps, and the value function again admits the form in terms of the scale function. We also conduct 
a series of numerical experiments using the spectrally negative Levy process with hyperexponential jumps so as 
to verify the optimality of the proposed strategies and also the efficiency of the proposed algorithm. We refer the 
reader to, among others, [1, 4, 5, 6, 15, 16, 22, 25] for optimal stopping problems of spectrally negative Levy 
processes. 

The rest of the paper is organized as follows. In Section 2, we review the spectrally negative Levy process and 
the scale function and then solve the one-stage problem. In Section 3, we extend it to the multiple-stage problem. 
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In Section 4, we verify the optimality and efficiency of the algorithm through a series of numerical experiments. 
Section 5 concludes the paper. 

2. One-stage Problem 

Let (O, J 7 , P) be a probability space hosting a spectrally negative Levy process X = {X t : t > 0} characterized 
uniquely by the Laplace exponent 

1 



(2.1) V(/3) := E u e^ Al = c/3 + -a 2 f3 2 + / (e""* - 1 + f3zl {0<z<1} ) U(dz), p G R, 

L J 2 J(0,oo) 

where c G M, c > and II is a Levy measure concentrated on (0, oo) such that 

(2.2) / (1 A z 2 )Ii(dz) < oo. 

Here and throughout the paper P^ is the conditional probability where Xq = i£l and E x is its expectation. We 
exclude the case X is a negative subordinator (decreasing a.s.) and we shall further assume that the Levy measure 
is atomless: 

Assumption 2.1. We assume that II does not have atoms. 

This section considers the one-stage optimal stopping problem of the form (1.2) where the supremum is taken 
over the set (or a subset) of stopping times with respect to the natural filtration F = {Ft)t>o generated by X. We 
assume the running payoff function / to be locally bounded and increasing. The stochastic process X models the 
state of the project and the monotonicity of / means that it yields higher benefits when X is high. Typically one 
assumes f(x) = e x — 5 as in (1.1) and this is clearly a special case of this model. Regarding the terminal reward 
function g, we consider two cases: (i) when g is a sum of exponential functions (Assumption 2.2 below) and (ii) 
when g is a general decreasing and concave function (Assumption 2.3 below). 

The results discussed in this section are applications of Egami and Yamazaki [17] and will be extended to the 
multiple-stage problem in the next section. Let S be the set of all [0, oo] -valued F-stopping times and define for 
any r G S, 



(2.3) u(x,T) = u(x,T;f,g) := E a 



r e- rt f(X t )dt + e- rT g(X T )l {T 
Jo 



<oo} 



x G 



After a brief review on the scale function and the results of [17], we shall solve under Assumption 2.2 below that 

u(x) := sup u(x, t). 
res 

We then obtain under Assumption 2.3 below a weaker version of optimality 

u(x) := sup u(x, r), 
res 

over the set of all first down-crossing times, 

S :={T A :Ae R}, 

where 

(2.4) t a := inf {t > : X t < A} , AG R, 
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with inf = oo by convention. This form of optimality is often used in real options and also in the field of 
corporate finance and credit risk as exemplified by Leland's endogenous default model [23, 24]. In practice, a 
strategy must be simple enough to implement and it is in many cases a reasonable assumption to focus on the set 
of stopping times of threshold type as in (2.4). Because S C S, it is clear that u >u. For the rest of the paper, let 
h±(x) := ±h(x) V 0, x 6 M, for any measurable function h : M. \- > M.. 

2.1. Review of scale functions and Egami and Yamazaki [17]. For any spectrally negative Levy process, there 
exists a function called the (r-)scale function 

W {r) r > 0, 

which is zero on (— oo, 0), continuous and strictly increasing on [0, oo), and is characterized by the Laplace trans- 
form: 



/ e - SX W^( X )dx = — , 8 > $ r , 

Jo Vis) - r 

where 

$ r := sup{A > : V(A) = r}, r> 0. 

Here, the Laplace exponent ip in (2.1) is known to be zero at the origin, convex on M + ; $ r is well-defined and is 
strictly positive whenever r > 0. We also define 

Z^'\x):=l + r[ W (r \y)dy, xGR, 
Jo 

which is also called the scale function. As we shall see below, the pair of scale functions and play 
significant roles in our problems. 

As in Lemmas 8.3 and 8.5 of Kyprianou [21], for each x > 0, the functions r i-> (x) and r h-> (x) can 
be analytically extended to r G C. Fix a > and define ip a {'), as the Laplace exponent of X under P a with the 
change of measure ^| = exp(aX t — ip(a)t), t > 0; as in page 213 of [21], for all f3 > —a, 

(2.5) MP) ■= (aa 2 + c - ^ u( e - au - l)n(du))/? + -a 2 (3 2 + J (e^ u - 1 + pul {0<u<l} ) e - au U(du). 

If W a and Z a are the scale functions associated with X under P a (or equivalently with ip a (-))- Then, by Lemma 
8.4 of [21], 

(2.6) Wj r -^ a »(x) = e- ax W^(x), x>0. 
In particular, by setting a = $ r (or equivalently r = tp(a)), we can define 

<f>r 

which satisfies 



(2.7) Wq> r {x) := w£\x) = e-^W«(i), x £ 



[•CO 1 

Jo ip(0 + $r)-r 



The smoothness and asymptotic behaviors around zero of the scale function are particular important in our 
analysis. We summarize these in the remark given immediately below. 
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Remark 2.1. (1) Assumption 2.1 guarantees that is C 1 on (0, oo). In particular, when a > 0, then 
is C 2 on (0 , oo). Fore more details on the smoothness of the scale function, see Chan et al. [12]. 
(2) As in Lemmas 4.3 and 4.4 of [22], 



0, unbounded variation 
- , bounded variation 



WM'(0+) :=lim^W'(x) 
x^o 



\, (J > 

OO, (7 = ant/ 11(0, oo) = oo 

r+n(o,oo) ^ compound Poisson 



ft- 

where fj, := c + Jj ^ z LT(dz), which is finite when X is of bounded variation. 

In [17], we have shown that a candidate optimal stopping time can be efficiently identified using the scale 
function. Define the expected payoff corresponding to the down-crossing time (2.4) by 

(2.8) ua(x) := u(x , ta) , i,i G M, 

which necessarily equals g(x) for x < A. By combining the compensation formula for Levy processes and the 
resolvent measure written in terms of the scale function, this can be written in a semi-explicit form. Let 



(2.9) 



(2.10) Q f (x;A): 
and 



poo 

* f (A) := / e-^yf(y + A)dy, 4el, 
Jo 

' J2w( r \x-y)f(y)dy, x > A, 



x < A, 



p^ A := LT(du) £ e-^ z (g(z + A - u) - g(A))dz 

rco ru+A 

(2.11) =/ n(d«) / e-^y- A \g(y-u)-g(A))dy, Ael, 

JO J A 

rco /•uA(x-A) 

v W(x) := U(du) W ( - r) (x-z-A)(g(z + A-u) - g(A))dz, x > A. 

Jo Jo 

These integrals are well-defined if 

/•CO 

(2.12) / e-*""|/(y)|dy < oo, 
Jo 

rco 

(2.13) g^C 2 , and / U(du) max \g(() - g(A)\ < oo, AgR. 

Jl A-u<C<A 

If these are satisfied, we can write ua{x) as in (2.8) for x > A as the sum of the following three terms: 



Z^(x-A) - ^W {r \x-A) 



ri(s;A) :=g(A) 

(2 " 14) r 2 (x;A) :=W^(x-A)p^ A -^ A (x), 

T 3 (x;A) :=W (r) {x-A)^ f {A)-e f {x;A) 
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Egami and Yamazaki [17] obtained the first-order condition that makes 8ua(x) / OA vanish and showed that it 
is equivalent to the continuous fit condition ua{A+) := linx^ ua(x) = g(A) when X is of bounded variation 
and to the smooth fit condition u' A {A+) := \\m. x ^Au' A {x) = g'(A) when a > 0. Recall that X is of bounded 
variation if and only if a = and 

(2.15) f (1 Az)II(dz) < oo. 

J(0,oo) 

It has been shown that 

(2.16) u A {A+) =g(A) + W ( - r \0)A(A), Ael, 
where 

2 

(2.17) A(A) = A(A; f,g) := ~g(A) - ^g'(A) + pfy + V f (A), AeK. 

In view of Remark 2.1 (2), for the unbounded variation case, continuous fit holds whatever the choice of A is, 
while, for the bounded variation case, it holds if and only if A(A) = 0. 

Furthermore, it has been shown, on condition that there exists some 5 > satisfying 

(2.18) / Il(dit) sup \g(A + £) -g(A + £-u)\ < oo, 
we have 

(2.19) A UA ( x ) = _ e *r(*-A) w ^ (a ._ A)A(A)) X>A} 

where W$ r is defined in (2.7). It is known that W$ r is increasing and hence, if A(A) is monotonically increasing, 
the down-crossing time ta for the A with A(A) = becomes a natural candidate for the optimal stopping time. 

2.2. Exponential Case. We first consider the case where g admits the form: 

N 

(2.20) g(x) = K~Y^ Cie a ' x , x e R, 

i=l 

for some constants K € K and Cj, a, > 0, 1 < i < N. We assume without loss of generality that Oj 7^ dj for 
z 7^ j. Because g is bounded on (—00, B) for any £> £ R, (2.13) and (2.18) are clearly satisfied. For /, we need a 
technical condition so that (2.12) is guaranteed. 

Assumption 2.2. Suppose 

(1) /(•) is continuous, increasing, /(— 00) := linx^-oo f(x) > —00 and J °° e~® rV f + (y)dy < 00; 

(2) <?(•) admits the form (2.20) for some K £ R awe/ strictly positive constants en and c%, 1 < z < iV, swc/i 
?/j<3? aj 7^ (ijfor any i 7^ j. 

Remark 2.2. For special case r > maxi<j<jv V'( a i) or equivalently <I> r > maxj = i ..jy dj, we can simplify the 
problem and rewrite it so that g is a constant and f is an increasing function. 
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With Assumption 2.2, we simplify (2.17) using 

By the monotonicity of ip, w r (a) > for any a > 0. The proof of the following lemma is given in Appendix A.l. 
Lemma 2.1. For every A G R, we /jave 

(2.22) A(A) = -—if + ^ cie^^tu^oi) + * f (A). 



In view of (2.22) above, the function A(A) is clearly continuous and increasing. Therefore, if lim^-oo A(A) < 
< lim^oo A(A), there exists a unique root A* G R such that A(^4*) = 0. Otherwise, let A* = — oo if 
lim j4 |„ 00 A(A) > and let A* = oo if lim^oo A (.A) < 0. 

Remark 2.3. (1) For any stopping time r G S, recall Assumption 2.2 {1) and noticed [JJ" e~ rt f(— oo)dt] = 
- E x [e -rr ]). Therefore, with f(x) := f(x) - /(-oo) > andg(x) := g(x) - 

u(x) = \-mpu(x,T;f,g). 

r res 

Because g is decreasing, there exists a unique A E [-co, oo] such that g(x) < <J=^ x G (A, oo). 
Because f is nonnegative, it is easy to see that it is never optimal to stop on (A, oo) and hence we can 
also write 

u[x) = h sup u(x,T\},g + ). 

r res 

(2) Applying (2.17) to the pair f and g, because / is uniformly positive, 

N 



A(A; lg) = ~L.(K- ^-^) + E Cie^Cai) + VP^A) 



AT „, s JV 



i=l j=1 



where the second to last inequality holds because, by the convexity of tp on [0, oo), w r {a,i) > If 
— oo < A < oo, because A is the unique root of — -^-g(A) = 0, we must have A* < A. If A = — oo, we 
have A(A) > -^-g(A) > for all A G M am/ Zierace A* = -oo. 
(3) Except for the case g is a constant, because g decreases to — oo, A < oo. Hence (2) implies A* < oo. 



When — oo < A* < oo, the corresponding expected value becomes (see, e.g., Exercise 8.7 (ii) of [21]) 

N 
i=i 

+ W {r \x- A*)V f (A*) -Q f (x;A*). 



u A *(x) = K(z( r \x-A*)-^W {r \x-A*)^-Y,Cie aiX (z^^ 

i=i 
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By how A* is chosen in (2.22) and by (2.6), it can be simplified to 

N 

(2.23) u A *(x) = KZ^(x-A*)-Y,Cie a ^ x zi r :-^ a ^(x-A*)-Q f (x;A*). 

i=i 

The verification of optimality requires the following smoothness properties, whose proofs are given in Appendix 
A.2. 

Lemma 2.2. Suppose —oo<A*<oo. 

(1) u A *(x) isC 1 onR\{A*}. 

(2) In particular, when X is of unbounded variation, ua* (%) is C 2 on M\{A*}. 

Herein, we add a remark concerning continuous/smooth fit. The following remark confirms the results in [17] 
and further verifies that smooth fit holds whenever X is of unbounded variation even when a = 0. This observation 
only requires the asymptotic behavior of the scale function near zero as in Remark 2.1 (2). 

Remark 2.4 (continuous/smooth fit). (1) Continuous fit holds (i.e. ua*(A*+) = g(A*)) because, by (2.23), 
ZW(0) = Zj~^ (o<)) (0) = 1 andlim xlA * Q f (x;A*) = 0. 
(2) In particular, when X is of unbounded variation, smooth fit holds (i.e. u' A * (^4*+) = g'(A*)) because 

N 

u' A *(x) = KrW {r) {x - A*) -^2ae a ^(r - ip(a i ))w£-M ai)) \x - A*) 

i=l 

N N 

-J2^ie a ' x zi r -^(x - A*) - @'f(x; A*) ^ - £ a iCi e a ' A * = g'{A*) 

1=1 i=l 

thanks to W^(0) = w£~^ (a<)) (0) = 0, zj~* (ai)) (0) = 1 and lim xiA & f (x; A) = 0; see also the proof 
of Lemma 2.2. 

We now state the main results of this subsection. The proof is given in Appendix A. 3. 

Proposition 2.1. (1) If — oo < A* < oo, the stopping time ta* := inf {t > : X t < A*} is optimal over S 
and the value function is u(x) = ua* {x) as in (2.23) for all x & M. 

(2) If A* = oo, immediate stopping is always optimal and u(x) = g(x)for any x£l 

(3) If A* = —oo, it is never optimal to stop (i.e. t* = oo is optimal), and the value function is given by 



(2.24) 




•r(*-v)-W<rXy-xj)f(y)dy, 



where &' r is the derivative of <3? r with respect to r. 

2.3. For a general concave and decreasing g. We now relax the assumption on g and consider a general concave 
and decreasing function g. We also drop the continuity assumption on /. 

Assumption 2.3. Suppose 

(1) /(•) is increasing such that /(— oo) := lim :r |_ 00 f(x) > —oo and J °° e~® rV f+(y)dy < oo; 

(2) </(•) is twice-differentiable, concave and monotonically decreasing such that (2.13) and (2.18) hold. 



CONTRACTION OPTIONS IN SPECTRALLY NEGATIVE LEVY MODELS 



9 



Under this assumption, we see that A(A) is increasing. Indeed, we have 



d_ 
dA 



r , (T , , s . ( r ) 



r 



2 /*oo /*i 

- ^V) + / n(d«) / 

z Jo Jo 



( g '(z + A- U )-g'(A))dz. 



This is positive because g'(A) < and g"(A) < for any AgK and because by the concavity and monotonicity 
z + A — u<A for z<u and hence the integral is nonnegative. On the other hand, ^f(A) is increasing because 
/ is. Therefore, we again define A* in the same way as the unique root of A(A) = (if it exists). The proof of the 
following result is given in Appendix A.4. 

Proposition 2.2. Suppose Assumption 2.3. 

(1) When — oo < A* < oo, then ta* is optimal over S and the value function is given by 

2 

(2.25) u{x) = u A *(x) = g(A*)Z^(x - A*) + W^(x - A*)^g'(A*) - ^ A ,{x) -Q f (x;A*), x > A*. 

For x < A*, we have u(x) = g(x). 

(2) If A* = oo, immediate stopping is always optimal and u{x) = g(x)for any x € M. 

(3) If A* = — oo, then t* = oo is optimal over S and (2.24) holds. 



3. Multiple-stage problem 

In this section, we extend to the scenario the firm can decrease its involvement in the project in multiple stages 
as defined in (1.3). As in the one-stage case, we consider two modes of optimality: 



M 

(3.1) U {M \x):= sup Vff 

(rW,..,T(«))65 Mm= i 

M 

(3.2) U {M) {x):= sup J^E 2 

(rW,...,r(M)) € S Mm =l 



,. T ( m ) 

/ e- rt F m (X t )dt + e- rT<m) g m (X T(m) )l {T(m)<oo} 

J T (m-l) 



T (m) 



(m) 



r (m-l) 



e- rt F m (X t )dt + e- r ^ g m (X r(m) )l 



{ T ( m )<oo} 



for all x E M where we define := for notational brevity and the supremum is, respectively, over the set of 
increasing sequences of M stopping times, 

S M ■= {t^ g S, 1 < m < M : < < t^}, 

and over the set of increasing sequences of M down-crossing times, 

S M ■= {r (m) = r Am e S, 1 < m < M : A x > A 2 ■ ■ ■ > A M }. 

Clearly, S M C S M and hence U {M ^ < U^ M \ 
We first consider the case g m admits the form 



N„ 



(3.3) 



gm 



(x) :=K m -Y^ c mi e a ™ x 1 < m < M, 



i=l 



for some constants K m G K and c m i,a m i > 0, 1 < i < N m , and show the optimality in the sense of (3.1) as 
an extension of Proposition 2.1. We then consider a more general case where g m is twice-differentiable, concave 
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and monotonically decreasing and show the optimality over Sm as an extension of Proposition 2.2. Regarding the 
running reward function F, define the differences: 

fm ■= F m - Fm+i, 1 < m < M, 

with Fm+i = 0- As is also assumed in [9], we consider the case f m is increasing. Using the notation as in (2.3), 
we can then write for all x £ M 

M 

(3-4) U^ M \x)= sup Vu(x,r( m );/ m , 5m ), 

(r(D,...,T(M))6% m=1 
M 

(3.5) U^ M \x)= sup ^n(x,r( m );/ m , 5m ). 

(TW,...,r(M))eS M m=l 

In summary, we assume Assumptions 3.1 and 3.2 below for (3.1) and (3.2), respectively. 

Assumption 3.1. For each 1 < m < M, we assume that f m and g m satisfy Assumption 2.2. 

Assumption 3.2. For each 1 < m < M, we assume that f m and g m satisfy Assumption 2.3. 

As is clear from the problem structure, simultaneous stoppings (i.e. T& = • • • = t^+i a.s. for some k and /) may 
be optimal in case it is not advantageous to stay in some intermediate stages (i.e. stages k + 1, . . . , k + I). For this 
reason, define, for any subinterval X = {minX, minX + 1, . . . , maxX} C {1, . . . , M}, 

(3.6) g x := and f x := F minX - F max x+i, 

iex 

and consider an auxiliary one-stage problem (1.2) with g = gx and / = fx- Notice that Assumption 3.1 (resp. 
Assumption 3.2) guarantees that fx and gx also satisfy Assumption 2.2 (Assumption 2.3) for any X. Hence Propo- 
sitions 2. 1 and 2.2 apply. 
Let 

(3.7) A m {A) :=A(A;f m ,g m ), A e R, 1 < m < M, 

as the function (2.17) for (f m ,g m ). Because p^ +h2 ,A - Phi, A + Phl,A and ( A ) = ^hA A ) + ^h 2 ( A ) for 

any measurable functions hi and /i2, we see that 

(3-8) A X (A) := A(A;f x ,gx) = a(A; ^ f m , £ g m ) = £ A m {A) 

is increasing and corresponds to the function (2. 17) for (fx, gx)- In particular, under Assumption 2.2, this reduces 
to 

r Nm 

Ai(A) = ( " ^~ K m + ^2c mi e a ™ A nj r (a rm ) + ^ fm (A)J. 

m£l i=l 

Now let A* x be the root of A X {A) = if it exists. If lim^oo Ax(A) < 0, we set A* x = 00; if lirn^.oo A X (A) > 0, 
we set A x = —00. For simplicity, let A* m := A*^ m y for any 1 < m < M. Also define 

u\(x) :=u(x,T A ;fx,gx), x,AgR. 
With these notations, the following is immediate by Propositions 2.1 and 2.2. 
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Corollary 3.1. Fix any X and and consider the problems: 

u x (x) := sup u(x,t; fx, gx) and u x {x) := sup u(x, r; fx, gx). 
re5 res 

Suppose Assumption 3.1. 

(1) If -oo < A\ < oo, then 

N m 

ux(x) = u\{x) = Y, {KmZ^(x - A x ) -2<W^^ W) (i - AD) - e fl (x;A x ), 

mgl i=l 

tmcf stopping time r^j := inf {i > : Xt < AJ} is optimal. 

(2) /f^4j = °°> u l( x ) = gx{ x )for any i£l w/?/z optimal stopping time t* = 0. 

(3) /^^4j = ~~ °°' ' f iJ never optimal to stop, and the value function is given by 

/oo 
(<!>i r e-^~y)-w^(y-x))fx(y)dy. 
-oo 

Suppose Assumption 3.2. 

(1) If -oo < < oo, 

= ^.(x) = gx(A x )Z^(x - A x ) + W^(x - Ax)^g'x(A x ) - ^ (*) - e /z (x; AJ), * > AJ, 

with the optimal stopping time r^j. For x < A x , we have u{x) = g(x). 

(2) If A x = oo, ux{x) = gx( x )for any i£l with optimal stopping time r* = 0. 

(3) If A x = — oo, then r* = oo is optimal over S and (3.9) holds. 

3.1. Two-stage problem. In order to gain intuition, we first consider the case with M = 2 and obtain (x) and 
lJ( 2 \x) under Assumptions 3.1 and 3.2, respectively. Following the procedures discussed above, A* m E [— oo, oo], 
or the root of A m (A) = 0, is well-defined for m = 1, 2. As a special case of (3.6), 

(3.10) f 2 = F 2 , h = F 1 -F 2 = F 1 - f 2 , and f {1>2} = F 1 = h + f 2 . 

We shall consider the cases (i) A\ > A* 2 and (ii) A\ < A 2 , separately. For (i), we shall show that (t4*,ta*) is 
optimal. For (ii), we shall show that simultaneous stoppings are optimal. We first consider the former. 

Proposition 3.1. If oo > A\ > A 2 > — oo, then {ta*,ta*) is optimal; the value function is given x) = 

Sm=i 2 u a*( x ) andU^(x) = J2m=i 2 (x) under Assumptions 3.1 and 3.2, respectively. In particular, 
under Assumption 3.1 and if oo > A\ > A 2 > — oo, 

N m 

(3.11) U^(x)= £ (K m Z^\x-A*J-Y,c mi e a - x Z^ a ^\x-Al 



m=l,2 i=l 



W^(x-y)F 1 (y)dy- / W^(x - y)F 2 (y)dy. 
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Proof. Suppose Assumption 3.1 holds. By relaxing the constraint that < t^ 2 \ we can obtain an upper bound: 

U {2 \x) <U (2) (x) := SU P u(x,T^;f m ,g m )= £ «< m >/ 



m=l,2 T(m '^ s m=l,2 



where the last equality holds by Corollary 3.1. On the other hand, because ta* < ta* a.s. (hence (ta* , ta* ) G £2) 
thanks to A^ > A|, we have E/( 2 )(x) > [T j (ic), as desired. The same result holds under Assumption 3.2 by 
relaxing the constraint that A\ > A2 and noticing that (ta* , ta*) G £2- 
For the second claim, because A* > A\ and by (3.10), 

/x px p 

\ W<r\x-y)f m (y)dy= W^(x-y)F 1 (y)dy + I (x - y)F 2 (y)dy, 

and hence (3.1 1) holds in view of Corollary 3.1. 

□ 

Now consider the case —00 < A* < A\ < 00. 

Lemma 3.1. Suppose —00 < A\ < A\ < 00. Under Assumption 3.1 (resp. Assumption 3.2), the first optimal 
stopping cannot occur on (A 2 , 00); namely ifr*^ is the optimal first stopping time in the sense of (3.1) ((3.2)), 
then X T «(i) G (—00, A?;] a.s. on {t*W < 00}. 

Proof. The result is immediate when A\ = 00 and hence we assume A\ < 00. 

Suppose Assumption 3.1 holds, and in order to derive a contradiction, we suppose there exists some x > A\ at 
which it is optimal to stop. Under this assumption, the value function must satisfy 

(3.12) C/ (2) (x) = gi(x) + sup u(x,T;f 2 ,g 2 ) = 9i(x) + u A 2 )(x). 

res 2 

We shall show that this is in fact smaller than u^) (x) + u A 2 }(x ), which is the value obtained by (t^, t^) = 
(t A *,t a *) eS 2 C S 2 . By (2.16) and (2.19), 

(3.13) Ai(A) < (>)0=> -q^a } ^) > (<)0 VA<x, 

(3. 14) u« ( A +) = 91 (A) + (0) Ai (A) . 
By (3. 13)-(3. 14) and because A\ < A\ < x and Ai is increasing, 

u { })(x) > limn { 4 1} (x) = gi (x) + W (r) (0)Ai(x) > gUx). 

Regarding the last inequality, for the unbounded variation case, it holds because W^ r \0) = by Remark 2.1 (2). 
For the bounded variation case, it also holds because (3.13) and x > A\ imply Ai(x) > 0. Therefore, we get, by 
(3.12), U^(x) < u A )(x) + u A 2 )(x) leading to a contradiction. Because x is arbitrary on (A 2 ,oo), we have the 
claim. The same contradiction can be derived under Assumption 3.2 because (ta* , ta* ) G S 2 . □ 

The following lemma suggests under — 00 < A\ < A\ < 00 that the optimal strategy is the simultaneous 
stopping corresponding to the threshold level A^ 2 y, which is the value that makes A{ lj2 } = Ai + A2 as in (3.8) 
vanish. 
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Proposition 3.2. Suppose — oo < A* < A 2 < oo. 



(1) We have A\ < A* {12} < A\. 

(2) It is optimal to stop simultaneously and the value function is given by (x) = u J, (x) under As- 



sumption 3.1 and 



x = u 



{1,2} 
1 {1,2} 



[x) under Assumption 3.2. 



Proof. (1) Because both Ai and A2 are increasing, Ar lj2 } is increasing as well. Because A\ < A 2 , Ai(A 2 ) > 
and hence Ar 1)2 }.(A 2 ) > 0. Similarly, A{ lj2 }(^4*) < 0. The increasing property of A| X 2 } now shows the claim. 

(2) Under Assumption 3.1, for any pair of stopping times (t^\t^) G 5 2 , because F 2 = / 2 as in (3.10) and by 
the strong Markov property of the Levy process X, 



E" l > e 

m=l,2 



_ r _(m) 



9m (m))l| T (m) <00 } 



T m ~ r (2) 

e- rt F 1 (X t )dt+ / e- rt F 2 (X t )dt 



r(l) 



e~ rT 3i(X t(1) )1 {t(1)<oo} + / e- n F x {X t )dt 

'0 



(3.15) 



< E a 



+ l {T(1)<00} e— W E 

r-r(l) 



r(T (2)_ T (l) ) 



52(^ T (2))l{ r (2) <00 } + 



(1) 



r (2)_ T (l) 



t(rr,rW;Fi,^ 1 + 



L {r( 1 )<oo} 
0h 



-(1) 



ffl (X T(1) )l |r(1)<oo} + / e-^F^Xjdt + l fT(1 



r(l) 



L {r( 1 )<oo| e " supu(A r( i),r;/ 2 ,5( 2 ] 
re-S 



Similarly, under Assumption 3.2, (3.15) also holds for any (t^\t^) £ 1S2 by replacing S with S in the 
expectation of the third term. This together with Lemma 3.1, shows (x) (resp. (x)) is less than or equal to 



(3.16) 



sup u(x,r;Fi,g 1 + u x A J 



{2K 

2 ' 



under Assumption 3.1 (Assumption 3.2) where S(A 2 ) is the set of r G 5 (r G 5) such that X T G (—00, A|] a.s. 

r 2 i 

on {r < 00}. Because, for any r G 5(^4 2 ), (A r ) = g2(X T ) a.s. on {r < 00} and because F\ = /{i i2 } as m 
(3.10), (3.16) equals 



-fl 2\ 

sup u(x,r;/{ 1>2 i,5r lj2} ) < supu(x, r; / {1;2} , #{i, 2 }) = «V ( x ) 

reS(A*) re-S t 1 ' 2 } 



by Corollary 3.1 and because S(A* 2 ) C S. Namely, f/( 2 )(x) < M^ 1 ,' 2 ^ (x) under Assumption 3.1 and U^ 2 \x) < 

{1'2} 

■fl 21 11 21 

u A i (x) under Assumption 3.2. These in fact hold with equality because?/^*' (x) is attained by (ta* {1 2} > TA^ 2} ) G 

□ 
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3.2. Multiple-stage problem. We now generalize to the multiple-stage problem and solve (3.1)-(3.2) or equiva- 
lently (3.4)-(3.5) with M > 3. For 1 < m < M, let 

M 

(3.17) UM{X):= sup Vn(x,r( fc );/ fc , 5fc ), 

( r (m) i ... iT CM) )e5M _ m+1 k=m 
M 

(3.18) U^\x):= sup ^u(x,r( fe );/ fc , 5fe ). 

(T(™),.,r( M ))6S M _ m+ i fc=m 

In particular, f/W = [/^ and f/M = u[ M \ and by Corollary 3.1 

(3.19) u9f\x) = supu(x,r; /m,9m) = uV?* (x) and Zjj* f) (x) = supu(x,r; /m,9m) = u { ^? } (x), 

res re § m 

under Assumptions 3.1 and 3.2, respectively. The expressions for U^f\ and U^f \ can a l so De obtained as in the 
two-stage case. 

Given 1 < m < M, let us partition {m, m + 1, . . . , M} to an L(m) number of (non-empty) disjoint sets 
X m := {l(k; m), 1 < k < L(m)} such that 

{m, m + 1, . . . , M} = X(l; m) U • • • U l(L(m); m) 

where, if L(m) = 1, 1(1; m) = {m, . . . , M} and, if L(m) > 2, 

Z(l; m) := {m, . . . , m, m - 1}, 

X(Z; m) := . . . , n z>m - 1}, 2 < Z < L(m) - 1, 

X(L(m); m) := {n L(m) _ 1]m , . . . , M}, 

for some integers m < ni >m < • • • < Ji2,( m )_i, m < M. We consider the strategy such that, if and Z are in the 
same set, then the k-th and Z-th stops occur simultaneously a.s. 

We shall show that (3.17) and (3.18), for any 1 < m < M, can be solved by a strategy with some partition 
X* m := {I*{k; m),l<k< L*(m)} satisfying 

Ar*(l;m) > " " > ^J*(L*(m);m)> 

where A* x is defined as in (3.8) for any set X. The corresponding expected value becomes 

L*(m) L*(m) 

(3.20) ^?l>):= E ^ T ^ (fc:m) ;/x*(M^z*(W = J] u ^r2) (x) ' 

fc=l k=l 

whose strategy is given by for any m < n < M, 

T *(n) _ ta*^ k } for the unique 1 < k < L*(m) such that n G X*(fc; m). 

We shall show that (3.20) is optimal, i.e. Um^ = U^O* under Assumption 3.1 and Um = U^X under 
Assumption 3.2 for any 1 < m < M. Moreover, X^ can be obtained inductively moving backwards starting from 
X* M with L*(M) = 1 and X*(l; M) = {M}. For the inductive step, the following algorithm outputs X^_ x from 
X^ for any 2 < m < M. By repeating this, we can obtain the partition X* ; the resulting U^-l as in (3.20) becomes 



the value function Z7 (A/) = U{ 



(M) 
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Algorithm T* m _ x = Update(Xj^ , m) 
Step 1: Set i = 1. 
Step 2: Set 




{m - 1} UX*(l;m) U • • ■Ul*(i - l;m), i > 2. 



Step 3: Compute At and 

(1) if % = L*(m) + 1, then stop and return 1*^ = {X*(l;m - 1)} with L*(m - 1) = 1 and 
X*(l;m - 1) = {m- 1, . . . , M}; 

(2) if A* ± > A* XHi . m) , then stop and return X* m _ x = {l*(k;m - 1),1 < k < L*(m - 1)} with 
L*(m - 1) = L*(m) - i + 2 and 

(3.21) X*(l;m-l)=X and X*(Z;m - 1) = X*(Z + i - 2;m), 2</<L*(m-l); 

(3) if A% < At tt .,, set i = i + 1 and go back to Step 2. 



The role of the algorithm is in words to extend from n(= M— m+l)-stage problem to n+l(= M— m+2)-stage 
problem. The idea is similar to what we discussed in the previous section how to extend from a one-stage problem 
to a two-stage problem. When a new initial stage is added, the corresponding threshold value At is first calculated. 
Depending on whether its value is higher than that of the subsequent stages or not, simultaneous stoppings may 
become optimal. For n larger than two, we must solve it recursively by keeping updating the set X, or the set of 
the first (simultaneous) stoppings, as given in this algorithm. If At. is low, the strategy of the new initial stage 
may naturally depend on the strategies of all the subsequent stages. Unlike the extension to the two-stage problem 
which only needs to take into account the strategy of the stage immediately next, it needs to reflect the strategies 
of all subsequent stages. 

We prove the following under Assumption 3.1 for the optimality (3.1). As is already clear after the detailed 
discussion on the two-stage case, only a slight modification is needed for (3.2) under Assumption 3.2. 

Lemma 3.2. In view of the algorithm above, suppose Assumption 3.1 and fix 2 < m < M. Given thatX^ satisfies, 
for every 1 < I < L* (m), 

L*(m) 
k=l 

and is used as an input in the algorithm. Then, we have the following. 
(1) At the end of Step 2,ifl<i< L*(m), 

L*(m) 

(3.23) U^_\(x)<sup u(x,r;f ± ,g ± ) + V u(x,^ T (A^ (k , m) ); fa*(k;m), 9i*(k;m)) 

res 1 " " r~. 

k=% 

where ~d p {A) := v + TA°O v for any v £ S and A £ M with the time-shift operator Of, andifi = L*(m) + 1 

TT {M) / \ / f x {m-l,...,M} / \ 

T £S {m-l,...,M} 
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(2) Letl^ n _ 1 be produced by the algorithm. For any 1 < I < L*(m — 1), 

L*(m-1) 



(3.24) 



Tj(M) 



minl*(i;»n— 1) > 



X* (fc;m — 1) 



Proof. (1) We shall show by mathematical induction. 

(Base-step) Suppose i = 1. By our assumption (3.22) and by an argument similar to (3.15), 



— TT 

e 9r, 



1 (X T )lt T<oo} + / e- rt F^i(X t )di + e^tfW(X T )l 



{r<oo} 



sup 

TS.S 



L*(m) 



t(x,r;/ m _i,g m _i) + ^ «(^,^T(^i»(fc ;m ));/z*(fe;m),3z*(fc;m)) 



fc=l 



Now (3.23) holds because X = {m — 1} for i = 1. 

(Inductive-step) Now we assume (3.23) fori = j < L*(m) — 1, i.e.,X = {m— 1}L)X*(1; m)U- • -UX*(j-l; m) 
and 

L*(m) 

(3.25) [/^(x) < sup u(x,T-J ± ,g ± ) + ^ tl(x, l?r(^x.(fe ;m )); fx*(k;m), 9X*{k;m)) , 



k=j 



and show that it will hold for i = j + 1. 

Because when At > v the algorithm stops at j and never returns to Step 2, we suppose here that 

A*~ < Aj,,.. m y In view of the right-hand side of (3.25), if there exists some x > &x*t-. m \ at which it is optimal to 



stop, then the value function becomes g x (x) + ^^x*(i-m) ^ our assumption (3.22). Using the same reasoning 
as in Lemma 3.1, this is in fact smaller than vP A * (x) + U^. n/ s (x). Hence is is never optimal to stop on 
{Ax*(j-m) } 00 ) ^ or tne optimization problem on the right-hand side of (3.25) (see also the proof of Proposition 3.2). 

Now let 5(^2* (j. m )) be the set of all stopping times at which X G (-co, ^x*(j-m)] aS- ^ or a ^ T G ^4«(]-m))' 
we have r = M^z*(j;m)) a.s. and hence u(x, ^r{A* x , {j . m) )] fx*(j;m), 9x*{j;m)) = u i x i r ; h*(j;m), 9x*(j-m))- 
Therefore (3.25) implies 



U^_\(x) < sup 



< sup 

res 



L*{m) 

U(X,T; f±LIX*{j;m)'9±uX*(j; m )) + U ^ Xl ^ T ( A X*(k;m))i fx*(k;m), 9x*(k;m)) 

k=j+l 

L*{m) 

u(x,T;f XuXHj . m) ,g XuXHj . m) ) + ^ U ( X iMAz* { k;m)yjx*(k;m),9l*(k;m)) 



k=j+l 

Hence, (3.23) holds for i = j + 1, as desired. This proves (1) by mathematical induction. 
(2) When the algorithm stops, it is either (i) i = L* (m) + 1 at Step 3 (1) or (ii) At > A 



X* (i;m) 



at Step 3 (2). 



(i) Suppose i = L* (m) + 1. In this case, X = {m — 1, ... , M} and, by (3.23), 



U 



mini* 



(l-m-l)( x ) = U m-l( x ) ^ su P^(^ r ;/{m-l....,Ai},5{m-l....,M}) 



u 



res 



{m-l,...,M} 
l* ( 

{m-l,...,Af} 



which in fact holds by equality because the right-hand side is attained by (ta* ,...,ta* ) € 

' M ^ fc ^ V /1 {m-l,...,M} ' ' A {m-1,...,M} J 



~>M-m+2 



CS 



M-m+2- 
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(ii) Suppose the algorithm exits at i with A*^ > A*x*u. m y By (3.23), we have 

L*(m) 



U^Jxix) < sup u{x,r;f ± ,g + sup ^ u(x, ^r{^x*{ k -,m))'^ fx*(k;m), 9x*{k,m))- 



Regarding the second supremum of the right-hand side, the strategy {t"'\ min X*(i, m) < I < M}, defined by 
= ^r(^4j»( fc . m )) for the unique i < k < L*(m) such that I G X*(k;m), is feasible (or in <SM+l-minX* (i,m)) 
for any stopping time r G 5 and therefore 

L*(m) 

SUp ^ "(^^rC^Z^fe^));/!*^),^^)) 



r&5 fc= 



< _ _ ^ sup Yl < x M k) ;fk,9k) = uW t x ). 



/- T (minX*(i,m)) rt*'))^,,,, . _»,. . , . \ 

Hence, we obtain abound U^f_\(x) < u\* (x) + ^i.^jW- This together with (3.21) and (3.22) shows 



L*(m) L*(m-1) 
fc=i fc=l 

This holds by equality because the left-hand side is attained by a feasible strategy defined by 1* n _ x . This shows 
(3.24) for case I = 1. 

On the other hand, for any 2 < I < L*(m - 1), by (3.21) and (3.22), 

L*{m) L*(m-1) 

U (M) (x) - U {M) (x) - V u x * ik ' m) (x) ~ V u 1 "^™- 1 ) fx) 

fc=i+i-2 fc=Z 

which guarantees (3.24), as desired. 

□ 

Using Lemma 3.2 as an inductive step, the main theorem is immediate. Indeed, (3.22) holds trivially for M by 
Corollary 3.1. By applying the algorithm M — 1 times, we can obtain (3.1) for M — 1, M — 2, . . . , 1. 

Theorem 3.1. Let {Z„; 1 < m < M} fee produced by the algorithm. 

(1) Under Assumption 3.1, for every 1 < m < M and 1 < i < L*(m), 

L*{m) 

AM) i s _ X*(k;m) 



TjyM) ( r ) _ v f T ) 



particular, 



(3.26) ^ M Hx) = ^ (1;1) (x)=E^! ) W- 



fe =l 



18 



K. YAMAZAKI 



(2) Under Assumption 3.2, for every 1 < m < M and 1 < i < L*(m), 



L* (m) 
k=i 



In particular, 

0.2D ^ (m) (-)=^ (1;1 )(-)=e^: ) w- 

k=l 

4. Hyperexponential jump diffusion case and numerical examples 

In this section, we consider spectrally negative Levy processes with i.i.d. hyperexponential jumps and provide 
numerical examples. If a Levy measure has a completely monotone density, it can be approximated arbitrarily 
closely by that of a compound Poisson process with hyperexponential jumps (see, e.g., [18, 19]). In a related work, 
Asmussen et al. [3] approximate the Levy density of the CGMY process by a hyperexponential density. Herein, 
we will use the explicit expression of the scale function obtained by [18]. 

4.1. Spectrally Negative Levy Processes with Hyperexponential Jumps. Let X be a spectrally negative Levy 
process of the form 



Aft 

(4.1) X t - X = p,t + oB t - ^ Z n , 0<t<cx 

n=l 



Here B = {B t ; t > 0} is a standard Brownian motion, N = {N t ;t > 0} is a Poisson process with arrival rate A, 
and Z = {Z n ; n = 1, 2, . . .} is an i.i.d. sequence of hyperexponential random variables with density function 

J 

K z ) ■= ^2PjVje~ VjZ , z>0, 

3=1 

for some pi H h pj = 1 and < r\\ < ■ ■ ■ < rjj < oo. Its Laplace exponent (2.1) is given by 

1 J s 
ip(s) = fis + -a 2 s 2 - X^pj , sGl. 

j=i + s 

In particular, this reduces to a Brownian motion when J = and to (spectrally negative) exponential jump diffusion 
when J = 1. 

For our examples, we assume a > 0; see [18] for the case a = 0. In this case, there are J+1 negative solutions 
to the equation tp(s) = r and their absolute values {Ci,r',i = 1, . . . , J + 1} satisfy the interlacing condition: 
< £i r < rji < ^2,r < ■ ■ ■ < VJ < £ J+i,r < oo. For this process, the scale functions are given by 



(4.2) 

J+1 J+1 r -, 

(x) = ^Ci [e* pa! - e~^x and (x) = 1 + r ^ Q — (e*^ - l) + — (e~^ x - 1 
i=l i=l L r ?v 
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for every x > where 



Ci : — As 



with A 



Uje{i,...,J} ( 1 



Hle{i,...,j+i}\{i} ( 1 



£i,7' 



1 < * < J + l. 



With $l c) := $ r - c and ^ := &, r + c, 1 < i < J + 1, for any c > 0, we have by (2.6) for any x > 



J+i 



e(<0„ 



(4.3) 



i=l 



J+l 



= 1 + (r - V(c)) E ( 



i=l 



1 , fl,( c )„ . 1 , A c )^. 



Thanks to their forms as sums of exponential functions, the value function we obtain below admits closed forms. 

4.2. Numerical results on the one-stage problem. We first consider the one-stage problem as studied in Section 
2. In our numerical examples, we consider two examples for g: 

(a) mixture of exponential functions: g( exp > as in (2.20); 

(b) linear function: g^ hn \x) := —ax + /3, x G M, for some a > and /3 G M. 

It is clear that g(' m ) satisfies Assumption 2.3 (2). Regarding /, we consider the following three examples: 

(i) simple function: f^ sim \y) := Z]_ 00<n<00 f^h n (y) f° r some constants • • • < /(" 2 ) < < /(°) < 
/(i) < /(2) < . . such that 

— oo < lim n |_ DO /W < lim^oo /W < oo and subdivisions I n := (Z n; Z n +i] 

of E; 

(ii) linear function with a lower bound: f^ hn \y) := b\[(y + 62) V 63] for some b\ > and 62, 63 G R; 

(iii) exponential function with an upper bound: f {exp) {y) := e( L ^ Ai3 for some L > and B 6i. 

These satisfy Assumption 2.3 (1) and in particular (ii) and (iii) satisfy Assumption 2.2 (1). Hence Proposition 2.2 
holds (or u = ua*) for any choice and in particular Proposition 2.1 holds (or u = u = ua*) for (a) with (ii) or (iii). 
In order to implement the optimal strategy, we first obtain A* using (2.17) (or (2.22) for (a)) and then compute 



the value function via (2.25) (or (2.23) for (a)). Toward this end, we need to compute p^ lin) A and ^' ( ' ra) A (x) 
and tyf(A) and @f(x; A) for (i)-(iii). The proofs of the following two lemmas are tedious but straightforward and 
hence omitted. 



„(r) 



Lemma 4.1. (1) For every A G M, we have 



(r) 
Pg(Un) >A 



A r 



5=1 



+ 



(2J For all x > A, 



(r) 



J+l J 

[x) = a\J2 C i^ZPj 

i=i j=i 



Vj(Vj+®r 



<5> r (x-A) _ -nj(x-A)\ 



11j(Vj ~ti,r) 
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Define, for any s < t and k > 0, 

J+l 

(k) , .X 

w 

1=1 



{k) (x-sVA) + _ <S> < r k) (x-tVA)+ 



This vanishes whenever [s,t] n [A, x] = 0. Also let for any s < t and a S M, 7^ a ^(s,t) := f^ye ay dy 



t _ l \ „at 
a a* l e 



Lemma 4.2. (1) For every A G M, 
1 



/(aim),--, $ 



— oo<n<oo 
6 3 



e -* r (Z n - J 4) + _ e ~$ r (i n+1 - J 4) + 



62 + -4 



1 _ e -$r[b 3 -(b2+A)]+ ) + 



-* r [6 3 -(62+A)] H 



+ e 



-*r[&3-(fe2+A)] H 



[h-(b 2 + A)] ±+ 1 



$ 2 



o-B „ P LA 



l-e-^'d-^ 



/« particular, for any sufficiently small 63, f j(k„) (A) = 61 
(2) For every x, A G K, 



fe 2 +A , _1 



e /(aim) (x;A) = / (n) ^°i(/n,Wi), 

— oo<n<oo 

e /(Kn ) (x; A) = 61 [fo 3 ^i°i( A > ( 5 3 - 62) V A A x) + &2wi°i((6 3 - 62) V A A x, x) 



J+i 



+ 61 Ci[e* pa, 7 ( ~* r) ((&3 - b 2 ) V A A x, x) - e -^ x ^\(b 3 - b 2 ) V A A x, z)], 



i=l 



8 /( exp) (z; A) = w®> A (A V (B/L) A x, x) + w£ A (A, A V A x). 

We are now ready to show our numerical results for the one-stage problem. Throughout the results given below, 
we use the common Levy process of the form (4.1) with \i = 1, a = 0.2, J = 2, p = [0.1, 0.9], rj = [1, 20] and 
A = 5. Namely, this is a hyperexponential Levy process with two kinds of negative jumps: frequent and small 
jumps and infrequent and large jumps. For g and /, we consider any combination of the following: 

(a) g = g ( ex P) with a = [0.1, 0.2, 0.3, 0.4] and c = [4,3,2,1]; 

(b) / = g( Hn ) with a = 1 and (3 = 0; 

and 

(i) / = ryfisim) with Ix = (_ QO) 0), I 2 = [0, 00), /« = -10 and /( 2 ) = 10; 

(ii) / = 7/(' m ) with b\ = 1 and b 2 = with sufficiently small ^3; 

(ifi) / = ryfiexp) W j th L = B = 1 . 

for the weight parameter 7 = 0, 0.05, 0.1. 

The results for (a) g = g( exp ^ and (b) g = g( /m ) are graphically shown in Figures 1 and 2, respectively. In each 
figure, we plot the function A(-) as in (2.17) and the value function u A * for each choice of /. As can be confirmed, 
the function A(-) is indeed monotonically increasing and hence the unique root A* of A(A) = can be obtained 



(i), 
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easily by the bisection method. Using these optimal threshold levels, the value functions are computed via (2.23) 
and (2.25) using Lemmas 4.1 and 4.2. 

We see that the value functions are differentiable even at the optimal threshold levels A* and this confirms the 
smooth fit as in Remark 2.4 because X is of unbounded variation with a > 0. In particular, for the case 7 = 
with no running reward (/ = 0), because the stopping value decreases to —00 as x increases, the value function 
converges to zero as x increases. For f( sim ^ with a discontinuity at zero, the value function turns out to have a kink 
at zero although it still appears to be differentiable. While these computations require bisection-type methods to 
obtain £j jr 's, $ r and A*, these can be solved instantaneously. 

In order to verify that these are indeed optimal, we focus on the case 7 = 0.05 and plot in Figure 3 the value 
function ua* in comparison with the expected values of "perturbed" strategies ua(-) for A = A* — 2, A* — 1, A* + 
1, A* + 2. Notice that these can be computed as the sum of T's as in (2. 14). For any choice of A, it is easy to see 
that ua is continuous as in Remark 2.4 but fails to be differentiable at A ^ A*. We can confirm in all six cases 
that ua* indeed dominates ua for A 7^ A* uniformly in x G M. This numerically verifies Propositions 2. 1 and 2.2. 

4.3. Numerical results on the multiple-stage problem. We now move onto the multiple-stage problem. We 
assume M = 3 for brevity and use for / and g the functions (a)-(b) and (i)-(iii) defined for the one-stage problem. 

We first verify (3.26) or the optimality over £3 for the case, for each 1 < m < 3, (a) g m = g( ex P) for some 
a m = (a m i)i<i<N m and c m = (c mi )i<i< Nm with N m = 4, and either (ii) f m = 7 m / (hn) with 61 = 1 and b 2 = 
with sufficiently small 63 or (iii) f m = j m f^ exp ^ with L = B = 1 for some 7 m > 0. 

We conduct a number of experiments for various values of {a m ,c m ,7 m , 1 < m < 3}. By using the al- 
gorithm given in Subsection 3.2, the optimal threshold levels A* = (A*^\ A*^ 2 \ yl*( 3 )) takes values among 
{A*, A* 2 , A* z , A*r x 2 ,, A*^ 2 3 ,,Arj 2 3 |} and satisfy one of the following four cases: 

Case 1: A*W = A*W = A*®; 
Case 2: A*^ > A*W = A*^; 
Case 3: A*^ = A< 2 ^ > A*^; 
Case 4: A*^ > A<^ > A*^ . 

Here we use a random random number generator to sample a m , c m and 7 m for each 1 < m < 3 and also 
the choice of f( lm )/fi ex p) until we attain each of Cases 1 to 4. The generated parameters and the corresponding 
threshold levels are summarized in Table 1 . In order to validate the optimality of the strategy (r^fi) , t a »(2) , t a *( 3 ) ), 

we compare in Figure 4 the value function with those of perturbed strategies (r^i) , t-y{t, , Tj(3) ), 1 < k < 6, where 

A k A k A k 

A k := (A*^\ A< 2 \ A*^) +5 k 

with 5 X := (1,0,0), 5 2 := (1,1,0),<5 3 J= (1,1,1),<5 4 := (0, 0, -1), 65 := (0,-1,-1) and 6 6 := (-1,-1,-1). 
It is clear that (tj(i) , r^-(2) , Tj(3) ) 6 1S3 C 1S3 because A^ > A^ > A^ by construction. Figure 4 suggests 
in all cases that the value obtained by (r A »(i) , t a *(2) , ^,(3)) dominates uniformly over x those obtained by the 
perturbed strategies. These results are indeed consistent with our main theoretical results as in (3.26). In view of 
Figure 4, we also observe that there are up to three kinks (at ^4*) in the value function although these still appear 
to be differentiable. This is again due to smooth fit as in Remark 2.4. The perturbed strategies on the other hand 
fail to be differentiable while they are still continuous. 
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A(-) for /(«"*») value function for 

FIGURE 1 . Plots of A and the value function ua* when g = g( ex P\ 
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-5 -4 -3 -2 -1 



X 

(a-i) g (ex Pl with /( sim ) 



7 



4 



3 



2 



1 








x 



(a-ii) with /^ in ) 




-5 -4 -3 -2 -1 

x 

(a-iii) with 



7 




1 1 1 1 1 ' 1 

-6 -5 -4 -3 -2 -1 

x 

(b-i) g( lin 1 with /( sim ) 




i 



-7 -6 -5 -4 -3 -2 

x 

(b-ii) g( hn 1 with /( Hn ) 

7 p 1 1 1 1 r 




X 



(b-iii) g^ with 



FIGURE 3. Verification of optimality: the value function ua* (solid black) and the stopping 
value g (solid blue) in comparison with ua,(-) for A = A* - 2, A* - 1, A* + 1, A* + 2 (dotted). 
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fl 


h 


h 


71 

/ J- 


72 


7a 
/ o 


Case 1 


lin 


exp 


lin 


0.2287 


0.0642 


0.7673 


Case 2 


lin 


exp 


exp 


0.2834 


0.8962 


0.8266 


Case 3 


lin 


exp 


lin 


0.8131 


0.3833 


0.6173 


Case 4 


exp 


exp 


exp 


0.4139 


0.3091 


0.2638 





ai 


02 


a3 


Case 1 
Case 2 
Case 3 
Case 4 


(0.34, 0.21,0.16, 0.43) 
(0.20, 0.42, 0.16, 0.44) 
(0.29, 0.12, 0.40, 0.27) 
(0.38, 0.40, 0.41, 0.25) 


(0.36,0.20, 0.41,0.25) 
(0.25, 0.30, 0.23, 0.36) 
(0.27, 0.23, 0.49, 0.04) 
(0.50, 0.10, 0.21,0.40) 


(0.32, 0.41,0.39, 0.32) 
(0.35, 0.29, 0.36, 0.01) 
(0.14, 0.11,0.01,0.40) 
(0.09, 0.50, 0.36, 0.18) 





Cl 


C2 


C3 


Case 1 
Case 2 
Case 3 
Case 4 


(4.75, 4.33, 4.99, 3.02) 
(3.37,0.59, 1.23, 2.73) 
(4.95,0.09, 2.67,3.13) 
(0.37, 0.97, 0.20, 2.79) 


(2.22,3.16, 1.12, 1.94) 
(2.19,4.07, 1.71, 2.81) 
(0.34, 3.42, 4.43, 0.69) 
(2.95,2.16,4.73,0.92) 


(0.30, 1.78,3.26, 0.71) 
(2.19, 1.62, 1.88, 1.98) 
(4.70,3.92,4.50, 1.09) 
(4.55, 3.75, 3.82, 2.49) 





A\ 


A* 


A* 
^3 


A* 


A* 
{2,3} 


A* 
{1,2,3} 




A *(2) 


A* (3) 


Case 1 
Case 2 
Case 3 
Case 4 


-3.37 
-0.75 
-2.57 
1.46 


-0.37 
-1.88 
-0.77 
-1.42 


-1.07 
-0.72 
-5.64 
-2.56 


-2.090 
-1.35 
-1.73 
-0.25 


-0.75 
-1.39 
-2.74 
-2.01 


-1.75 
-1.17 
-2.67 
-1.13 


-1.75 
-0.75 
-1.73 
1.46 


-1.75 
-1.39 
-1.73 
-1.42 


-1.75 
-1.39 
-5.64 
-2.56 


TABLE 1 . 


Parameters anc 


threshold levels 


for (3.26). 



We now move onto verifying (3.27) or the optimality over £3 for a wider class of / and g. We assume (b) 
9m = <7^ W ) for some a = a m and f3 = for each 1 < m < 3, (i) /1 = 7 1 / <s " ri ) with I\ = (— 00, 0), 
h = [0, 00), /W = -10 and /( 2 ) = 10, (ii) f 2 = 7 2 / (fcra) for h = 1 and b 2 = with sufficiently small 63 and 
( m ) /3 = 73/^ e2:p ^ for L = B = 1. We simulate the values of {a m ,j m , 1 < m < 3} until we obtain each of 
Cases 1 to 4 as described above. The generated parameters and threshold levels are summarized in Table 2 and the 
comparison with perturbed strategies are given in Figure 5. We again observe in all cases that the value obtained by 
(t a „ w , T A *( 2 ) , t a *(3) ) indeed dominates uniformly over x those obtained by the perturbed strategies. This verifies 
numerically (3.27) and we also observe the same smoothness properties as in Figure 4. 

5. Concluding Remarks 

In this paper, we studied a wide class of optimal stopping problems for a general spectrally negative Levy process 
and extended them to multiple- stopping. Our framework is applicable to a wide range of settings particularly in 
real option problems where the firm withdraw from a project in stages. Our analytical results suggest that the 
optimal solutions can be characterized by the threshold levels that make certain monotone functions vanish and the 
corresponding value functions can be expressed in terms of the scale function. Our numerical experiments suggest, 



26 



K. YAMAZAKI 




Case 3 Case 4 



FIGURE 4. Verification of optimality (3.26) for the multiple-stage case: the value function 
(solid) in comparison with the perturbed strategies. 

for the hyperexponential jump case, that these can be solved instantaneously with high precision. These tools we 
developed in this paper are highly valuable and can be used flexibly for analysis in real options and other fields of 
finance and industrial applications. 

There are several directions for future research. First, our results can be pursued for a general Levy process with 
both positive and negative jumps. While it makes the problem less tractable, it is expected that these can be done at 
least for the cases with analytical forms of Wiener-Hopf factors such as double exponential jump diffusion [20] and 
phase-type Levy processes [2]. Second, by using hyperexponential fitting, one can approximate any Levy process 
with completely monotone jumps (such as CGMY and variance gamma processes) by those with hyperexponential 
jumps as in Section 4. By calibrating with real financial and industrial data as in [3], one can conduct detailed 
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Oil 


Ct2 






I 6 


<Va 

/ o 


Case 1 


0.0596 


0.2687 


0.9867 


0.5824 


0.4210 


0.0921 


Case 2 


0.9119 


0.0149 


0.1567 


0.7722 


0.4754 


0.6809 


Case 3 


0.1874 


0.5317 


0.3550 


0.6003 


0.0850 


0.9224 


Case 4 


0.7862 


0.7231 


0.2788 


0.6535 


0.4897 


0.9729 





A\ 


A* 


A* 


A* 


A* 

^{2,3} 


A* 

^{1,2,3} 




A *(2) 


^4*(3) 


Case 1 


-4.53 


-4.7 1 


-1.56 


-4.61 


-2.84 


-3.47 


- 3.47 


- 3.47 


-3.47 


Case 2 


-3.32 


-6.69 


-6.20 


-3.93 


-6.52 


-4.15 


-3.32 


-6.52 


-6.52 


Case 3 


-4.17 


-2.12 


-4.87 


-3.47 


-3.42 


-3.79 


-3.47 


-3.47 


-4.87 


Case 4 


-3.31 


-3.58 


-5.60 


-3.41 


-4.14 


-3.70 


-3.31 


-3.58 


-5.60 



Table 2. Parameters and threshold levels for (3.27). 



empirical analyses on optimal stopping strategies and the value functions. Finally, it is an interesting extension to 
"swing option type" multiple- stopping with refracting time as in [10, 11] where there have to be some intervals 
between any two stoppings. 



Appendix A. Proofs 



A.l. Proof of Lemma 2.1. Substituting (2.20) in (2.11), 



(A.1) 



_(r) 
PgA 



N poo pu 

V Ci e a!A / LT(du) / e -*r*( e ai(*-u) - i)dz. 
i=1 Jo Jo 



(Case 1) First suppose ai ^ <P r for all 1 < i < N. Simple algebra gives 



(A.2) 
where 



HA) 



r 



N ( a i) 

-i=i 



Af W := r + -^-$ r + / n(d«) 



o 



1 - e 

By the definition of ip and <3? r , we rewrite as 

r 2 



e- au (l - e 



-(<J> r -a)u> 



<3? r - a 



aa 2 f°° 
r + — $ r + / n(dn 
* Jo 



a E K\{$ r }. 



(1 _ e -*.« _ $ r nl {ue(0il)} ) - e- au (l - e -(*r-«)«)_^_ + <P r ul {ue(0>1)} 

^ n(e- au - l)II(du)) $ r + y$r($r + a) 

n(d U )e- a " (l - e -^-a)u + (<1 > r _ a)«l {u6(0 , i; 



$ r - a 



i/j a (<f> r - a), 
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7 - 
6.5 - 



S 

a 4.5 ■ 



4 - 

3.5 - 

3 - 

2.5 — 



A +(1,0,0) 
A* +(1,1,0) 
A* +(1,1,1) 
A* -(0,0,1) 
A* -(0,1,1) 
A* -(1,1,1) 
■A* 



A +(1,0,0) 
A* +(1,1,0) 
A* +(1,1,1) 
A* - (0,0,1) 
A* - (0,1,1) 
A* -(1,1,1) 
A* 




Case 1 



Case 2 




A +(1,0,0) 
A* +(1,1,0) 
A* +(1,1,1) 
A* - (0,0,1) 
A* - (0,1,1) 
A* -(1,1,1) 
A* 



-4 -3.5 -3 -2.5 -2 




Case 3 



Case 4 



FIGURE 5. Verification of optimality (3.27) for the multiple-stage case: the value function U^> 
(solid) in comparison with the perturbed strategies. 



where the last equality holds by (2.5). On the other hand, tp a {^r — a) = ip{& r ) — VK a ) = r — ip(a); see page 213 
of [21]. Hence 

<S> r w r {a) = $r (r - tp(a)) = r ip a (^r ~ a), 
<P r — a <P r — a 

which shows for the case a» ^ <3? r for all 1 < i < N. 

(Case 2) Suppose aj = <£ r for some 1 < j < N (with ^ a,j for i ^ j by assumption). Take a sequence 
of (strictly) increasing sequence a - ' j aj = & r . Then a modification of (2.17) with aj replaced with a - is by 
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Case 1 



r 



AM(A) 

By the definition of w r as in (2.21), we have 



K+ Yl ^ A w T { ai ) + A Wr (aY ] ) + * f (A). 

l<i<N,i^j 



N 



lim A( m \A) = - —K + Y Ci e a * A zu r ( ai ) + V f (A). 



m\co 



i=l 



On the other hand, in view of (A.l), its integrand is monotone in a. Hence by the monotone convergence theorem 
and because g and g' are continuous in a, lim^oo A ( - m \A) = A(^4), and the proof is complete for Case 2. 

A.2. Proof of Lemma 2.2. Because g(x) is infinitely differentiable, the results are clear for x G (— oo,A*). 
Hence we show for x € (A*, oo). Because W^ T \y) is differentiable on y > as in Remark 2.1 (1), KZ^ r \x — 
A*) - Y^=i Cie a - X Zt'^ 11 ^ {x - A*) is twice differentiable. 

Regarding Qf(x; A*), as in the proof of Lemma 4.5 of [16], integration by parts, thanks to the continuity of /, 
gives 

r f W (r \x-y)f(y)dy 



f(y)(Z^(x)-Z^(x-y)) 



y=x 



y=A* 



f'(y)(ZW(x)-Z<r\x-y))dy 
f(x)(Z^(x) - 1) - f(A*)(ZW(x) - zW(x - A*)) - (f(x) - f(A*))Z^{x) + / ,/'(,/)Z" '(,- - ,/)«!,, 
-fix) + f(A*)Z^(x - A*) + f f{y)Z {r \x - y)dy. 



Hence 



@ f (x;A*) 



-f{x) + S{A*)Z<r\x - A*) + f f'iy)Z^ix - y)dy 

J A* 



Because Z^ r \0) = 1, 6/(x; A*) is differentiable with 



Q' f ix;A*) = fiA*)W^ix-A*)+ \ f(y)W^(x-y)dy. 



A* 



When X is of unbounded variation, because W^(0) = as in Remark 2.1 (2), 0/(x; ^4*) is twice-differentiable 
with 

&}ix; A*) = S(A*)WW'(x - A*) + |* S'{y)W^\x - y)dy. 

A3. Proof of Proposition 2.1. Suppose — oo < A* < oo. By directly using the results of [17] (Lemma 3.7 and 
Proposition 3.4), we obtain 

(C-r)uA*(x) + S(x) = 0, xe04*,oo), 
UA*ix) > gix), 
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where C is the infinitesimal generator of X applied to a sufficiently smooth function h, i.e., 

1 



Ch(x) = ch'(x) + -a h"(x) + / [h(x — z) — h(x) + h' (x)zlf 0<z< x\]H(dz). 

2 ./l! 

We shall show (£ — r)uA* (a;) + f(x) < on x G (— oo, ^4*) and verify the optimality. By the definition of ip, 



loo 

ca + -a a + 



Jo 



I + azl {0<z<1} ] U(dz) 



e ax *p(a) 



N 



i=l 



for any a > and hence we have 
(A.3) (£ - r)g(x) + f(x) = -rK + Ci e a ^(r - V(a,)) + f(x). 

By how j4* is chosen, 

(A.4) = -rK + ^2c i e a * A *<Z> r w r (a l ) + * r */(A*). 

Because / is increasing and x < ^4* 



iV 



i=i 



/(x)dy = /(x). 



r 

(A.5) $ r * f (A*)>$ r 

Jo 

It is also easy to see that for any 1 < i < N 

(A.6) e a * A * '$ r w r (ai) > e a * x (r - ^K)). 

Indeed, for the case r — ip(a%) > 0, we must have <3? r — <2j > and hence (A.6) holds by A* > x; for the 
case r — ip(a,i) < 0, the left-hand side is positive while the right-hand side is negative in (A.6); for the case 
r — ip(cii) = 0, the left-hand side is positive because ^>'($ r ) is, while the right-hand side is zero. Hence, by 
(A.4)-(A.6), (£ - r)u A *(x) + f(x) < holds. 

This result also holds for the case A* = oo. In this case, > -rK + Y^iLi Cie aiA <f> r m r (a,i) + ^ r ^f(A), for 
any A G M. Therefore (C — r)g(x) + f(x) < holds by the same reasoning as in (1) by simply replacing A* with 
A. 

Finally, we verify the optimality of ua* (x). Thanks to Lemma 2.2 and the continuous/smooth fit condition as in 
Remark 2.4, a version of Meyer-Ito's formula as in Theorem IV.71 of [28] (see also Theorem 2.1 of [27]) implies 
{e~ rt UA* (Xt) — Jq e~ rs (C — r)uA* (^s)l{jf a ^A*}ds, t > 0} is a local martingale. Namely, there exists some lo- 
calizing sequence {a m , m > 1} such that {e~ ritAam) u A *(X tAcTrn ) - f^*™ e- rs (C-r)u A *{X s )l {Xs ^ A *}ds, t > 
0} is a martingale for any m > 1. By the optional sampling theorem, for any stopping time r G S, we have for 
any L > 



E a 



-r(rAcr m AL) / y \ 



rAa m AL 



e rs (C - r)u A *(X s )l{x s ^ A *}ds =u A *(x). 



Because (C — r)uA* + / < 0, 



-r(rAa m AL) (y \ , 



rAo m AL 



e- rs f(X s )ds <u A *(x). 
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Moreover, because ua*(x) > g{x) and ua*(x) > lim^-oa ua(x) > Jo°° e r */( — °°)dt =: B uniformly by 
(2. 19), we have 



- r{TAamAL) (g(x TAamAL )vB) + 



-Aa m AL 



e- rs f(X s )ds <u A .(x). 



Because g V B £ [B, K V B], by the dominated convergence theorem 



lim lim E x 

mtoo L\oo 



r{ ™ L \g{X th(TmhL )y B)] =r[e-( S (I T )V5)l > E^ \e~ rT g{X T ) 



On the other hand, because /(— oo) > — oo and by the monotone convergence theorem, 

»TA(T m AL 



lim lim E 1 

m^oo Lfoo L Jq 



= lim lim E x 

m\oo L\oo 

This means E a 



e- rs f(X s )ds 

TAo m AL 

e- rs f + (X s )ds 



o 



lim lim E x 



TAtlmAL 



e~ rs f.(X s )d S 



= E x 


[/ 







e - r 7(x s )ds 



e rT g(X T ) + fj~ e rs /(X s )ds < u^. (x), for any arbitrary r G 5, as desired. 



It is now left to show for the case A* = — oo. Because 8ua{x) / dA < for any A € M as in (2. 19), there exists 



U-oo(x) := lim ua(x) = lim E x 
A^—oo Al—oo 



T~A 



e- n f(X t )dt 



o 



By Assumption 2.2 (1) and by the monotone convergence theorem, we can write 

l-T A 

E a 



u-oo(x) = lim ( E x 
A],— oo 



e- rt /+(^)dt 

" /"OO 

E* / e- rt /+(^)dt 
■/ o 









-rt 


/ 6 




.JO 




-E z 


u 




Jo 



= E x 


[I 







poo 

/ e- rt f(X t )dt 
Jo 



Now by Corollary 8.9 of [21], u-oo attains the proposed form. Because 8ua(x) / dA < for any A £ M, clearly 
-u_oo(x) > g(x) for any x£l. Moreover, because u_oo is attains by r* = oo, we have the claim. 

A.4. Proof of Proposition 2.2. (1,2) We first suppose x < A*. Then by definition ua*(x) = g(x). Because 
ua(x) = g(x) for any A > x, it is sufficient to show ua{x) < g(x) for A < x. This is indeed so because by 
(2.16), (2.19) and A(x) < due to x < A*, 

ua(x) < u x (x+) = g(x) + iy (r) (0)A(x) < g(x) = ua*{x), A < x < A*. 

This proves (2). For (1), suppose x > A*. Then by (2.19) ua*(x) > ua(x) for any A < x. For A > x, 
u A {x) = g(x) and by (2.16), (2.19) and A(x) > due to x > A*, 

u A *(x) > u x (x+) = g{x) + iy (r) (0)A(x) > g(x) = u A (x), A > x > A*. 

Therefore ua* (x) > ua(x) uniformly in x G M, as desired. 

The corresponding value function (for (1)) can be expressed as the sum of (2. 14): 

u(x) = g{A*)Z<r\x - A*) + W<r\x - A*) (~g(A*) + p% + *,(A*)) - <p^(x) - @ f (x; A*). 
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From the definition of A* that makes (2.17) vanish, 

u{x) = g(A*)zW(x - A*) + W<r\x - A*)°^g'(A*) - tp^x) - Q f (x; A*), 

as desired. 

(3) the proof is the same as that of Proposition 2. 1 (3). 
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